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ABSTRACT 



Context. The observation of rapidly variable very high energy (VHE) gamma-rays from non-aligned active galactic nuclei (AGNs), 
as reported from M87, proves challenging for conventional theoretical acceleration and emission models. 

Aims. Motivated by recent work on pulsar-type particle acceleration in M87 (Neronov & Aharonian 2007), we re-examine the cen- 
trifugal acceleration of particles by rotating jet magnetospheres in the vicinity of accreting supermassive black hole systems and 
analyze the energy constraints imposed for highly underluminous systems. 

Methods. The maximum Lorentz factor for centrifugally accelerated electrons in the presence of inverse Compton losses, and the 
associated characteristic variability time scale, are determined. Applications are presented for conditions expected to be present in the 
radio galaxy M87, assuming accretion onto the central black hole to occur in an advection-dominated (ADAF) mode. 
Results. We show that for a highly underluminous source like M87, centrifugally accelerated electrons may reach Lorentz factors up 
to y ~ (10 7 - 10 s ), allowing inverse Compton (Thomson) upscattering of sub-mm disk photons to the TeV regime. Upscattering of 
Comptonized disk photons results in a flat TeV spectrum L v oc v~" c with spectral index a c ^ 1.2. The characteristic variability time 
scale is of the order r^/c, which in the case of M87 corresponds to = 1.7 d for a typical light cylinder radius of tl - 5 r\. 
Conclusions. Centrifugal acceleration could provide a natural explanation for the challenging VHE emission features in M87. Our 
results suggest that some advection-dominated accreting (non-blazar) AGNs could well be observable VHE emitting sources. 

Key words, galaxies: active - galaxies: jets - radiation mechanism: nonthermal - gamma rays: theory - individual: M87 



1. Introduction 

The rapidly varying VHE gamma-ray flux, on time scales of days 
or less, observed in several AGNs implies a very compact VHE 
emission region of at most R < cAtd = 2.6 x 10 15 (A?/1 d) 5 
cm, where 6 is the bulk Doppler factor of the VHE emitting 
region. For blazar sources with their jets pointing almost di- 
rectly towards us (i.e., 6 ~ Ft, ~ 15), VHE variability on time 
scales of several hours may, in principle, be successfully ac- 
counted for by internal shock acceleration of electrons and in- 
verse Compton upscattering of soft photons (e.g., Mastichiadis 
& Kirk 2002). Yet, the fastest observed VHE blazar variations, 
on time scales of minutes (Gaidos et al. 1996; Aharonian et 
al. 2007; Albert et al. 2007), are generally difficult to under- 
stand within standard jet VHE emission models and are likely 
to require non-standard geometrical set-ups for their explana- 
tions (e.g., Salvati et al. 1998; Rieger 2004). In the case of non- 
aligned AGNs with 5 ~ 1, little can be gained from jet boost- 
ing, and it remains to be shown whether conventional, single ho- 
mogeneous SSC models may be flexible enough to reproduce 
the VHE characteristics, including variability on time scales 
Af 5 2 day, as observed in the radio galaxy M87 (Aharonian et 
al. 2006). Here we explore the possibility that centrifugal accel- 
eration of electrons, occurring in the vicinity of a sub-Eddington 
accreting black hole system, could via inverse Compton pro- 
cesses lead to variable VHE gamma-ray emission. Centrifugal 
acceleration of plasma flows by rotating magnetospheres has 
been widely discussed, both in the context of pulsar emission 
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models (e.g., Gold 1969; Machabeli & Rogava 1994; Chedia 
et al. 1996; Gangadhara 1996; Bogovalov 1997; Contopolous 
et al. 1999; Machabeli et al. 2005; Thomas & Gangadhara 
2007) and in the context of relativistic jet formation (e.g., 
Blandford & Payne 1982; Fendt 1997; Camenzind 1999; Meier 
et al. 2001). Detailed MHD simulations, for example, indicate 
that centrifugally-driven outflows from AGNs could reach bulk 
Lorentz factors ~ 10. It has been realized for almost a decade 
that the generated MHD jet topology could also allow efficient 
centrifugal acceleration of relativistic charged test particles to 
very high energies in sub-Eddington accreting black hole sys- 
tems (Gangadhara & Lesch 1997; Rieger & Mannheim 2000, 
henceforth RM00; Xu 2002; Osmanov et al. 2007). Here we ex- 
amine this issue in more depth for parameters relevant to the 
radio galaxy M87. 

2. Centrifugal acceleration of test particles 

2.1. Particle energization and acceleration efficiency 

We consider an idealized two-dimensional model topology 
where the magnetic field is assumed to rotate rigidly with a frac- 
tion of the rotational velocity of the black hole (e.g., Fendt 1997), 
and where the electric field component parallel to the mag- 
netic field line is screened by the magnetospheric jet plasma. 
A charged test particle, injected at the base and corotating with 
the field line (bead-on-wire motion), will then experience the 
centrifugal force and gain rotational energy while moving out- 
ward along the field (Machabeli & Rogava 1994; Machabeli et 
al. 1996; Gangadhara & Lesch 1997; RM00). The radial particle 
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motion can be most conveniently analyzed in the framework of 
Hamiltonian dynamics (RMOO) by noting that the Hamiltonian 
H for a particle (rest mass mo) on a relativistically rotating wire 
is a constant of motion, and is given by H = ymo c 2 (l-Q 2 r 2 /c 2 ), 
where y = 1/(1 - Q 2 r 2 /c 2 - r 2 /c 2 )^ 2 is the Lorentz factor, 
Q = c/tl the angular velocity of the field line and the light 
cylinder radius. The resultant equation of motion 



d 2 r dr dy 2 



(1) 



can be solved analytically yielding (Machabeli & Rogava 1994; 
RMOO) 



rit) = r L cn(A -Q.t) 



(2) 



for the time-dependence of the radial coordinate r, assuming a 
particle to be injected at ro with Lorentz factor yo, where cn is 
the Jacobian elliptic cosine and Aq is a Legendre elliptic integral 
of the first kind, and 



r(K0) 



l 



m(l -Q 2 r 2 (f)/c 2 ) 



(3) 



when expressed in terms of r(t), with m - I /(yo[l-Q. 2 r^/c 2 ]) 2 < 
1. If we neglect, for a moment, radiative losses and the break- 
down of the bead-on-the-wire approximation, a particle would 
reach the light cylinder r^ within a time t v = Ao/Q. < r^/c, 
where, as a consequence of the reversal of the centrifugal ac- 
celeration, it would change direction and move inward again 
(Machabeli & Rogava 1994). Knowing the dependence of y 
on r, we can easily determine the local acceleration time scale 
(RMOO) 



t - y -~ 

lace — . — 



c Vl - Cl 2 r 2 /c 2 



2D. 2 r Vl -m[l -Cl 2 r 2 /c 2 ] ' 



(4) 



which for large y approaches t^. = l/(2f2m'^ 4 y 1 ' 2 ) — 
16.6m~ 1/4 (10 6 /r) 1/2 ('-L/10 15 cm) [s]. The (local) acceleration 
time scale thus decreases with increasing y. Equation (01 implies 
that the particle Lorentz factor increases dramatically as a parti- 
cle approaches the light cylinder. Therefore, in many cases the 
characteristic linear size of our region of interest is much smaller 
than the light cylinder and the curvature radius of the field, a fact 
that may qualify the presumed straight field line approach. 

2.2. Efficiency constraints and maximum particle energy 

In reality, unlimited growth will be prohibited by radiative en- 
ergy losses, the breakdown of the bead-on-the-wire approxi- 
mation or the bending of the field line with increasing inertia 
(Gangadhara & Lesch 1997; RMOO). Inverse Compton upscat- 
tering of accretion disk photons, for example, leads to an en- 
ergy loss of the particle characterized by a time scale f coo i oc 1/y 
that decreases faster than f acc and so introduces a natural lim- 
itation. Indeed, only for highly underluminous AGN sources 
will centrifugal acceleration be sufficiently efficient to accelerate 
electrons well beyond Lorentz factors of one hundred (RMOO; 
Xu 2002). Assuming the inverse Compton scattering process to 
be approximately describable by the Thomson limit for a quasi- 
isotropic photon distribution with energy density t/ p h [erg/cm 3 ], 
we have t coo \ - 3 • 10 7 /[y t/ p h] s. Balancing acceleration by cool- 
ing (f acc = f coo i) for electrons thus gives a maximum electron 
Lorentz factor (see also Osmanov et al. 2007) 
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provided that the corotation condition can be satisfied for such 
a range of Lorentz factors. The latter qualification seems impor- 
tant, as it may well happen that the particle motion becomes so 
perturbed by radiation recoil that the bead-on-the-wire approxi- 
mation is no longer a useful concept. For highly underluminous 
AGN sources with, e.g. I/ ph ~ 0.01 erg/cm 3 and r L ^ 5 x 10 15 
cm (see below, for which Eq. (f5]l would otherwise imply 
Tmax ~ L3 X 10 9 m'/ 2 , this may in fact be the case. Indeed, from 
the bead-on-wire requirement that the radiation reaction force 
that results from inverse Compton scattering, i.e., F rat | P\c/c 
(with Pic the single particle Thomson power) should be (much) 
smaller than the Lorentz force Fl, one finds that achievable elec- 
tron Lorentz factors should be (much) smaller than 
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where fi(rjj is the magnetic field strength at the light cylinder 
radius. Yet, even for cases where radiative losses might be ne- 
glected, the breakdown of the bead-on-the-wire approximation 
(roughly occurring when the Coriolis force exceeds the Lorentz 
force) will prevent a particle from achieving infinite energies 
(RMOO). In the simplest case, this restricts achievable particle 
energies to Lorentz factors below Q 
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which implies lower Lorentz factors for protons than for elec- 
trons. In any case, once the Lorentz factors become too high, the 
inertia of the particles overcomes the tension in the field line, 
so that the line is swept back opposite to the sense of rotation, 
slowing down acceleration and introducing curvature radiative 
losses, thus ultimately preventing infinite energy growth. In what 
follows, it is assumed that achievable Lorentz factors y always 
satisfy the relation y < minlyj^ 



1 / max' I it 



3. Application to the radio galaxy M87 

3. 1 . Phenomenological background 

The nearby (distance ~ 16 Mpc) giant elliptical galaxy M87 
hosts one of the most massive black holes M BH - 3 x 10 9 M o 
(e.g., Marconi et al. 1997), with Schwarzschild radius r s = 
2GM BH /c 2 * 8.9 x 10 14 cm, and a prominent one-sided (kpc- 
scale) jet visible from radio to X-ray wavelengths (e.g., Marshall 
et al. 2002; see Ly et al. 2007 for a possible radio counter-jet de- 
tection). HST observations have revealed superluminal motion 
of jet components at ~ 0.5 kpc from the central black hole, in- 
dicative of bulk flow Lorentz factors Tf, ~ 6 and a jet orienta- 
tion of 6 ~ 19° to the line of sight (Biretta et al. 1999; but see 
also Ly et al. 2007 for larger radio 6), suggesting that M87 is a 
non-blazar jet source, characterized by only moderate Doppler 
factors. Superluminal radio features have also been detected in 
HST-1 located at around 100 pc (Cheung et al. 2007), although 
no superluminal motion has been found on small scales (Kovalev 
et al. 2007). HESS observations have recently shown M87 to be 



1 We are grateful to Osmanov et al. (2007) for making us aware of 
an incorrect conclusion in RMOO. Although formula (19) in RMOO is 
correct, an inconsistent set of parameters has been used to estimate yj^ x 
for the applications presented. This becomes relevant for highly under- 
luminous sources, where achievable Lorentz factors can be much higher 
than previously concluded in RMOO. 
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a rapidly variable (observed time scale of ~ 2 days) TeV emit- 
ting source, yet with a relatively low (isotropic) TeV luminosity 
of several times 10 40 erg/s (Aharonian et al. 2006). The total 
nuclear (disk and jet) bolometric luminosity of M87 has been 
estimated to be of order L\, \ ~ 10 42 erg/s or less (Reynolds et al. 
1996; Owen et al. 2000), indicating that M87 is a highly under- 
luminous source with l e < 3 x 10~ 6 , where l e = Lboi/^Edd and 
Lfidd is the Eddington luminosity. This has led to the proposal 
that M87 is a prototype galaxy, where accretion occurs in a two- 
temperature, advective-dominated (ADAF) mode characterized 
by an intrinsically low radiative efficiency (Reynolds et al. 1996; 
Camenzind 1999; Di Matteo et al. 2003). Magnetic flux dragged 
inwards may then build up a rigidly rotating, dipolar magneto- 
sphere, along which disk plasma can be centrifugally accelerated 
to (bulk) outflow Lorentz factors of Ft, - (5-10). The generated 
light cylinder scale is likely to be of order ~ 5 r s (Camenzind 
& Krockenberger 1992; Fendt 1997; Camenzind 1999; cf. also 
Fendt & Memola 2001 for higher if differential rotation is 
important). The magnetic field lines in global (quasi force-free) 
MHD wind solutions are radial near the black hole horizon, but 
asymptotically collimated into a cylindrical structure beyond the 
light cylinder, typically on radial scales ~ 10 This seems con- 
sistent with high frequency VLBI observations in M87, indicat- 
ing a jet that forms with opening angles ~ 60° at the jet base 
(Ly et al. 2007), as expected in MHD models, and a jet radius of 
~ 50 r s close to the origin (Rrichbaum et al. 2006). It has been 
proposed recently that efficient pulsar-type particle acceleration 
may occur in such an environment (Neronov & Aharonian 2007, 
henceforth NA07). 



3.2. Implications for particle acceleration 

Let us consider the implications of these findings for the cen- 
trifugal acceleration of particles in M87: 

(1) Firstly, even for the most limiting case where all of the 
observed bolometric luminosity is assumed to originate within a 
region tl, so that the energy density of the radiation field may 
be approximated by U p h = Lb i/(4nr^c), Eq. © would imply 
that Lorentz factors yj^ ax ^ 10 7 m^ 2 can be achieved, allowing 
Thomson upscattering of infrared (> 0.01 eV) photons to the 
TeV regime. In reality, this case is certainly over-restrictive, as it 
assumes that neither the observed jet nor the disk regions beyond 
tl make a significant contribution to the bolometric luminosity 
output, which we consider unlikely. Indeed, if the relevant lu- 
minosity is an order of magnitude or more smaller, as expected 
in the ADAF scenario (see estimate Lr below), the maximum 
Lorentz factors implied by Eq. © will be at least two orders of 
magnitude higher. This suggests that electron Lorentz factors up 
to y ~ (10 7 - 10 s ) may be well possible (cf. Eqs. JH-gl). 

(2) If accretion in M87 indeed occurs in an ADAF mode, 
the emitted disk spectrum will range from the radio up to the X- 
ray regime and beyond: the radio part is produced by synchro- 
cyclotron emission of thermal electrons (T e ^ 5 x 10 9 K), the 
optical/UV/X-rays arise via inverse Compton scattering of ra- 
dio soft photons, and the hard X-rays are due to bremsstrahlung 
and multiple Compton scattering (Mahadevan 1997; Narayan et 
al. 1998; Yi 1999). The ADAF equipartition magnetic field for 
M87 is of order B eq ~ 2.5 x 10 4 m 1/2 (r/r. 5 ) _5/4 G, where m is 
the accretion rate in units of the Eddington rate (cf. Yi 1999). 
For the inferred Bondi accretion rate m = /«/, ~ 1.6 x 10~ 3 (Di 
Matteo et al. 2003) this gives B(r s ) ~ 10 3 G, suggesting possible 
(radial) jet magnetic field strengths close to the light cylinder of 
Bifh) ~ (10 - 50) G. The highest radio emission in an ADAF 



is generally produced in the innermost region of the accretion 
flow. For M87, the characteristic synchrotron (peak) frequency 
becomes 



v 5 (r)-4xl0 13 m!, /2 (-) f T ' „ ) Hz, 
V ; * \r s ) \5 x 10 9 K/ 

and the associated radio luminosity Lr ~ v s Ll is given by 

x 8/5 , J x 21/5 / y ^7/5 



(8) 



/s, (9) 



h \10 3 / \5xl0 9 K7 \10 n Hz/ s 

where (xm/10 3 ) ~ 1 denotes the dimensionless synchrotron self- 
absorption frequency (Yi & Boughn 1998). On the light cylinder 
scale (rL - 5r s ) this implies a peak frequency v s {r]_) =2 2 x 10 11 
Hz and a luminosity L R ~ 2 x 10 40 erg/s. Thomson upscattering 
(v ~ y 2 v s ) of these mm soft photons by centrifugally acceler- 
ated electrons with y up to ~ (10 7 - 10 8 ) will thus, in princi- 
ple, lead to VHE photons with energies up to ~ (0.1 - 10) TeV. 
Comptonization of cyclosynchrotron soft photons adds a power 
law tail to the disk spectrum above v s , i.e., L v ^ L Vt {v/v s )~ ac 
with power index a c = -lnr es /lnA (Mahadevan 1997). For 
M87 with ihb = 1.6 x 10~ 3 , viscosity parameter a = 0.3 and 
T e =5xl0 9 K, we obtain an electron scattering depth r es ^ 0.04 
and an amplification factor A 15.7, so that a c ^ 1.2. 

(3) Suppose that during an active state, test particles are in- 
jected with 70 ~ 2 at a constant rate Q and accelerated up to a 
threshold y\, < y™ x < y^ x , above which they are considered to 
leave the centrifugal acceleration process due to the breakdown 
of corotation. The differential particle energy distribution n(y) 
along a field line would satisfy the simplified transport equation 



dn d 
~dt + dy 



y 
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^cool 

and U, 



"I 5(y h -y) = Q6(y- y ) , (10) 

^"esc 



with f coo i oc and ? acc oc y~ x l 2 [cf. Eq. |@)]. Above injection, 
the steady-state distribution in the acceleration region thus be- 
comes 



n{y) oc 1 — 



ic 

/max 



y 



-3/2 



H{y b -y), 



(11) 



i.e., a power law distribution with index -3 /2 for y^ <K y l ^ x . The 
emergent (singly scattered, Thomson) inverse Compton spec- 
trum ./'ic(v) from such a hard electron distribution would fol- 
low a power law jic(v) oc y~ ff with index a = 0.25 for v <k 
4y^v,. Integrating Eq. (fTT~b over y gives the number of parti- 
cles along a field line N — Q/O.. Electrons, escaping quasi- 
monoenergetically with y& ~ 10 7 from the acceleration mech- 
anism and encountering the Comptonized disk photons (v > v s ), 
can Thomson upscatter them to the TeV regime, producing a 
power law-like energy distribution above 4-y^v s with index a c ~ 
1.2, consistent with the value 1.22 ± 0.15 derived for the HESS 
2005 observations of M87 (Aharonian et al. 2006). 

(4) The number of escaping particles per unit time is n e ~ 
"(Ti)7fc/ T esc- Thus, within some time smaller than the cooling 
time (At = pt coo \,P < 1), we accumulate Nb ~ n e At particles that 
can IC upscatter Comptonized disk photons. We can roughly es- 
timate the associated TeV luminosity from L\c ~ Nb Pic, where 
Pic = l3cr T cylUph is the single particle Compton power per 
unit volume. This gives 



- M, '(^f(p)(S"' 



S/s. 



(12) 



To achieve a Compton luminosity comparable to the observed 
TeV luminosity of Lxev — 3 X 10 40 erg/s (Aharonian et al. 2006), 
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we thus need N ~ 3x 10 jp particles along field lines. Denoting 
the relevant acceleration volume by AV ~ 77 r£ Ar, with charac- 
teristic length scale Ar = \y/(dy/dr)\ ~ (yo/yb) r L and 77 < 1, 
the corresponding kinetic energy density n{ji,)y 2 h m e c 2 /AV (for 

prj ~ 10~ 5 ) is still well below the energy density B 2 /(Sn) of 
the magnetic field, suggesting that the presumed (quasi force- 
free) MHD field structure is still a valid approximation (cf. also 
Osmanov et al. 2007). 

(5) In principle, TeV gamma-rays can be strongly attenuated 
due to photon-photon pair production in the background disk 
photon field. The narrow dependence of the cross-section cr yr 
on the product of photon energies implies that VHE photons of 
energy E interact most efficiently with infrared background pho- 
tons of energy ei R (ITeV/E) eV. The optical depth t for a y- 
ray photon in a background field of infrared luminosity L 1R and 
size Ri R thus becomes (cf. NA07) 

indicating that due to its low bolometric luminosity M87 could 
be well transparent to VHE gamma-rays, even if almost all of 
the observed infrared luminosity Lir 10 41 erg/s (Whysong & 
Antonucci 2004) is (somewhat unrealistically) taken to be pro- 
duced on a scale Rir ~ r\ ~ 60 r s . Note that even if t would 
become larger than one, y-rays from the last transparent layer 
are still able to escape, so that the VHE flux would not simply 
decrease exponentially by exp(-r), but only by a factor of ~ r 
(NA07). 

(6) The number of electrons escaping quasi- 
monoenergetically from the acceleration mechanism in the 
vicinity of the light cylinder is of order Nh- Once these 
energetic particles encounter non-vanishing perpendicu- 
lar and/or the turbulent plasma magnetic fields, they can 
produce synchrotron emission arising as L Y oc y 1 / 3 be- 
low, and decaying exponentially above the peak frequency 
v syn ~ 50(y h /5 x 10 7 ) 2 (B sinor/lG) MeV with a total lumi- 
nosity of order L syn ~ P syn Nh ~ 0.06 Lie (B sin a) 2 , where P syD 
is the single particle synchrotron power. In order to satisfy the 
restrictions imposed by the existing (yet non-contemporaneous) 
upper limit on the M87 flux in the EGRET energy band above 
100 MeV (e.g., Reimer et al. 2003), the effectively encountered 
fields should be smaller than ~ 1 Gauss. This seems consistent 
with independent estimates suggesting a strength of the random 
field component close to the black hole of below one Gauss 
(NA07). The overall spectral energy distribution in M87 is then 
likely to consist of a number of different contributions, involving 
also other leptonic (Georganopoulos et al. 2003; NA07) and 
perhaps even hadronic (Reimer et al. 2004) processes. If so, then 
no straighforward X-ray-TeV correlation might be expected. 

(7) As shown above, accelerating particles up to the light 
cylinder typically takes a time r^/c, suggesting a characteristic 
variability time scale for M87 of t v ^ ~ ^ 1 .7 days, well 
consistent with the observed TeV time scale of Af ~ 2 days, a 
fact that may further validate the assumptions of the presented 
model. 

4. Conclusions 

VHE radiation from low-luminous, non-blazar AGN jet sources 
like M87 could provide an ideal test laboratory for the analysis of 
particle acceleration processes close to the supermassive black 
hole event horizon. In blazars with their relativistic jets point- 
ing towards us, most of these traces are likely to be masked by 



strong relativistic beaming effects, while for luminous quasars 
internal absorption of gamma-rays becomes dominant. Based 
on a simple toy model we have shown that efficient centrifu- 
gal acceleration of electrons in the vicinity of the light cylinder 
could provide a natural explanation for variable (time scale of 
one day) VHE emission with a hard inverse Compton spectrum 
as observed in M87. Our models fits well with other evidence 
for advection-dominated accretion in M87 and may indeed be 
regarded as providing some further corroboration for the pres- 
ence of such modes in highly underluminous AGNs. 
As always, there are a number of subtleties whose impact on the 
presented results need to be explored in more details including 
general relativistic effects, anisotropic scattering modifications, 
quasi rigid rotation and plasma instabilities. The extent to which 
our conclusions might be affected may require fully relativistic 
modelling. Yet, given the demonstrated potential of centrifugal 
acceleration and our current understanding of relativistic jet for- 
mation, this may represent a program worth pursuing. 

Acknowledgements. Discussion with and comments by John Kirk, Christian 
Fendt and Karl Mannheim are gratefully acknowledged. 
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